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Abstract. Sudden changes in stratospheric chlorine species 
in the polar northern atmosphere, caused by the Solar Proton 
Events (SPEs) of 17 and 20 January 2005, have been investi- 
gated and compared with version 4 of the Whole Atmosphere 
Community Climate Model (WACCM4). We used Aura Mi- 
crowave Limb Sounder (MLS) measurements to monitor the 
variability of CIO, HC1, HOC1 and Michelson Interferometer 
for Passive Atmospheric Sounder (MIPAS) on ENVISAT to 
retrieve CIONO 2 . SPE-induced chlorine activation has been 
identified. HC1 decrease occurred at nearly all the investi- 
gated altitudes (i.e., 10-0.5 hPa) with the strongest decrease 
(of about 0.25 ppbv) on 21 January. HOC1 was found to 
be the main active chlorine species under nighttime condi- 
tions (with increases of more than 0.2 ppbv) whereas both 
HOC1 and CIO enhancements (about 0. 1 ppbv) have been ob- 
served at the polar night terminator. Further, small CIO de- 
creases (of less than 0.1 ppbv) and CIONO 2 enhancements 
(about 0.2 ppbv) have been observed at higher latitudes (i.e., 
at nighttime) roughly above 2 hPa. 

While WACCM4 reproduces most of the SPE-induced 
variability in the chlorine species fairly well, in some particu- 
lar regions discrepancies between the modeled and measured 
temporal evolution of the abundances of chlorine species 
were found. HOC1 changes are modelled very well with re- 
spect to both magnitude and geographic distribution. CIO 
decreases are reproduced at high latitudes, whereas CIO 
enhancements in the terminator region are underestimated 
and attributed to background variations. WACCM4 also re- 


produces the HC1 depletion in the mesosphere but it does 
not show the observed decrease below about 2 hPa. Finally, 
WACCM4 simulations indicate that the observed CIONO 2 
increase is dominated by background variability, although 
SPE-induced production might contribute by 0. 1 ppbv. 


1 Introduction 

Since the early seventies it is known that energetic solar 
protons are able to affect the neutral chemical components 

of the polar atmosphere, such as [NO x ] = [NO] + [NO 2 ] 

and [HO x ] = [H] + [OH] + [HO 2 ] via ion chemical reactions 

(e.g., Swider and Keneshea, 1973; Crutzen et ah, 1975; 

Solomon et ah, 1981). These phenomena are important be- 
cause the NO x and HO x enhancements lead to efficient cat- 

alytic cycles of ozone destruction. The short chemical life- 
time of mesospheric hydroxyl (OH) radicals causes a rela- 
tively brief atmospheric impact from the increased HO x dur- 
ing Solar Proton Events (SPEs) (e.g., Damiani et ah, 2010a). 
A short term ozone depletion occurs from the enhanced SPE- 
produced HO x , followed by a fast recovery. On the other 
hand, the longer chemical lifetime of NO x , in particular dur- 
ing the polar winters (when photolysis does not occur), can 
cause noticeable ozone depletion in the upper stratosphere 
lasting up to several months (e.g., Randall et ah, 2005; Jack- 
man et ah, 2009). For this reason, previous efforts have been 
mainly focussed on this important relationship (i.e., NO x vs. 
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O 3 ). Recently, thanks to new satellite observations, many 
studies have dealt with the impact of SPEs (e.g., Jackman 
et al., 2005; Rohen et al., 2005; Semeniuk et al., 2005; 
Egorova et al., 2011), in particular the so-called Halloween 
SPE of October-November 2003, on other chemical com- 
ponents (e.g., HNO 3 , N 2 O 5 , CIONO 2 ) (e.g., Lopez-Puertas 
et al., 2005a,b; Funke et al., 201 1). The variability of chlorine 
species is a further interesting consequence of solar protons 
entering the Earth’s atmosphere. This has been investigated 
only in a few studies (i.e., Solomon and Crutzen, 1981; von 
Clarmann et al., 2005; Winkler et al., 2009) but none of these 
reported on variations in all the main components (i.e., HC1, 
CIONO 2 , CIO and HOC1) of the chlorine family during a sin- 
gle SPE. This is an essential condition to correctly analyse 
SPE effects. 

Hydrogen chloride (HC1) and chlorine monoxide (CIO) 
are the main contributors to the Cl y family ([Cl y ]=[Cl] + 
[CIO] + [HC1] + [CIONO 2 ] + [2C1 2 0 2 ] + [HOC1] + [OCIO] + 
[BrCl] + [2CL]). HC1 is the most important stratospheric 
reservoir species for Cl. CIO carries most of the reactive 
chlorine in the stratosphere during daytime. Their monitoring 
is very important since both gases can be considered a signif- 
icant indicator of the effectiveness of the Montreal protocol 
agreement (e.g., Froidevaux et al., 2006; Jones et al., 2011). 
On the other hand, a significant interplay between NO x and 
C10 x is provided by CIONO 2 , which is another important 
reservoir species. Also hypochlorous acid (HOC1) has re- 
markable importance for its catalytic cycle of mid-latitude 
lower stratospheric ozone loss (Lary, 1997), and because it 
acts as a link between chlorine and HO x chemistry. How- 
ever, very few HOC1 observations from satellites have been 
available until recently. The reason for this is twofold. First, 
this is due to the fact that HOC1 is not a very useful indi- 
cator of the atmospheric chlorine content from human made 
CFCs; second, there is also a technical difficulty due to the 
weak nature of HOC1 transitions in the microwave/millimeter 
wave region and in the infrared windows, and the low con- 
centration of HOC1 and low signal/noise ratio. 

In the upper stratosphere, the chlorine catalytic cycles 
of ozone destruction have the greatest efficiency around 
40 km (Lary, 1997), where SPE-induced changes of chlo- 
rine species are expected. There, they can contribute to the 
short-term O 3 depletion occurring after a SPE (von Clarmann 
et al., 2005), which is principally caused by enhanced NO x 
and HO x . In this regard. Salmi et al. (201 1) have found that, 
although not connected to SPEs, the ozone decrease of late 
February 2009 occurring at polar northern latitudes between 
about 30 and 50 km, was due to activation of halogen chem- 
istry. Although SPE-induced changes of chlorine species are 
expected to have an impact only on the short-timescale SPE- 
induced O 3 variability, the study of chlorine changes could 
be useful to validate chemical schemes in current models. 

The influence of SPEs on the chlorine family is in- 
direct and it occurs via the increase of HO x and NO x . 
SPE-induced enhancement in HO x is the main factor 
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leading to chlorine changes. Indeed, in prior studies 
(e.g., von Clarmann et al., 2005; Damiani et al., 2009) chlo- 
rine variations, recorded under SPE conditions, were used as 
a proxy of HO x changes. Moreover SPE-induced enhance- 
ment in NO x is important with respect to CIONO 2 . The pio- 
neering work of Solomon and Crutzen (1981) proposed sev- 
eral chemical mechanisms to explain the impact of the SPE 
of August 1972 on polar atmospheric composition. Two re- 
actions included in those schemes relevant to this study are: 

C10 + N0 2 + M-> CIONO 2 +M (Rl) 

and 

HC1 + OH Cl + H 2 0. (R2) 

The first experimental confirmation (by satellite data) of 
the importance of Reaction (Rl) under SPE conditions was 
reported by Lopez-Puertas et al. (2005a) and von Clarmann 
et al. (2005). The former examined the CIONO 2 in the 
context of NO y components (i.e., HNO 3 , N 2 O 5 , CIONO 2 
and NO x ), whereas the latter in the context of chlorine 
species (i.e., CIO, HOC1, CIONO 2 ). In both works, data from 
Michelson Interferometer for Passive Atmospheric Sounder 
(MIPAS), recorded during the 2003 Halloween storm, were 
used. In particular, von Clarmann et al. (2005) showed an 
enhancement of CIO and HOC1 immediately after the SPE, 
which was interpreted as being due to the fast reaction: 

CI + O 3 -> CIO + O 2 (R3) 

(with Cl arising from Reaction R2), and a subsequent HOC1 
formation by: 

CIO + H0 2 -» HOC1 + 0 2 . (R4) 

In this context a key role is played by HO 2 . The SPE-induced 
HO 2 enhancement may lead to a shift of the balance be- 
tween CIO and HOC1 towards the latter, so, particularly un- 
der nighttime conditions, a decrease of CIO could also occur 
(Funke et al., 201 1). 

There is still a further important pathway which can sig- 
nificantly reduce the ClO/HCl ratio in the upper stratosphere 
and, hence, the sensitivity of ozone to chlorine: 

CIO + OH -» HC1 + 0 2 . (R5) 

Indeed, although the major products of the OH + CIO reac- 
tion are the active HO 2 and Cl radicals, even a small branch- 
ing ratio may lead to less ozone depletion by chlorine species 
since Reaction (R5), faster than Reaction (R2), converts an 
active form (CIO) into a stable reservoir species (HC1) (e.g., 
Lipson et al., 1997). 

Some preliminary confirmation of the relevant role of Re- 
action (R2) has been shown by Damiani et al. (2009) using 
Microwave Limb Sounder (MLS) observations recorded dur- 
ing the SPEs of January 2005. They show an HOC1 strato- 
spheric enhancement accompanied by a similar decrease in 
HC1 that clearly demonstrates SPE-induced chlorine activa- 
tion. 
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The conversion of mesospheric HC1 into active chlorine 
was extensively examined by Winkler et al. (2009) using 
HALOE data recorded between 62° and 69° N during the 
SPE of July 2000. In this prior study, the temporal evolution 
of SPE-induced changes was partially reported since the em- 
ployed dataset included only three days after the event. Their 
simulations showed that the magnitude of the HC1 loss can- 
not be explained solely by the gas phase Reaction (R2) and 
model runs with parameterized production of NO x and HO x 
underestimate the HC1 loss. Indeed several negative ions re- 
act with HC1 and form negative chlorine species (e.g.. Cl - ) 
mainly in the mesosphere. Then, they can react with H to re- 
form HC1, or with positive ions releasing, e.g.. Cl, CIO. The 
increasing ion production rates during a SPE can change the 
concentration of charged chlorine species and it results in a 
final net HC1 depletion (additional to the depletion induced 
by OH) and an increase of other neutral chlorine species (e.g. 
Cl, CIO) (see Winkler et al., 2009). Adding to the model re- 
actions of HC1 with negative ions improves the agreement 
between simulation results and observed data in the meso- 
sphere (above about 55 km). However, a consistent model- 
data discrepancy in the stratosphere (where the impact of the 
negative chlorine chemistry is very small) remains (see also 
Winkler et al.,2011). 

The present paper deals with chlorine changes occurring 
in the Northern Hemisphere (NH) polar region during Jan- 
uary 2005 SPEs. Changes in the concentrations of O 3 , HO x , 
NO x , HNO 3 and in temperature caused by the SPEs of 17 
and 20 January 2005 have been reported in recent studies 
(e.g., Seppala et al., 2006; Verronen et al., 2006, 2007, 2011a; 
Klekociuk et al., 2007; Damiani et al., 2008; Seppala et al., 
2008; Orsolini et al., 2009; Becker and von Savigny, 2010; 
Jackman et al., 2011). In this paper we study the changes that 
occurred in a larger set of chlorine species (i.e., HC1, CIO, 
HOC1 and CIONO 2 ) by using MLS and MIPAS data and 
Whole Atmosphere Community Climate Model (WACCM) 
version 4 simulations. In our approach we did not consider 
ion chemistry (as in Winkler et al., 2009), so we indirectly 
tested its importance in the lower mesosphere. On the other 
hand, we studied the temporal evolution of chlorine changes 
by a larger dataset from polar orbit satellites, that provides a 
more representative picture of the polar cap region with re- 
spect to the prior work (i.e., Winkler et al., 2009). 

Recent modelling results (Jackman et al., 2008), compared 
with MIPAS data recorded during the Halloween storm, re- 
vealed good agreement with HOC1 data but poor results 
for CIO. These difficulties in reproducing SPE-induced CIO 
changes have been confirmed in Funke et al. (2011). More- 
over, model-data discrepancies are also evident when repro- 
ducing SPE-induced changes in CIONO 2 (Jackman et al., 
2008; Funke et al., 2011). A limitation in these prior studies 
(i.e., Jackman et al., 2008; Funke et al., 2011) has been the 
lack of measurements of the most important chlorine reser- 
voir (i.e., HC1). On the other hand, when HC1 has been inves- 
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tigated (Winkler et al., 2009), satellite observations of other 
components were missing. 

Some difficulties in studying and reproducing the effects 
of the SPEs of October-November 2003 (i.e., von Clarmann 
et al., 2005; Jackman et al., 2008; Funke et al., 2011) were 
likely due to the early winter vortex being partially illumi- 
nated. Since this paper deals with chlorine changes occurring 
in the NH during January, when the vortex is well established 
and coincident with the terminator, better modelling results 
are expected. 

Some evidence has been provided (Funke et al., 201 1) that 
the model-data discrepancy in CIO could arise also from the 
underestimation of its “background” (i.e., CIO VMRs be- 
fore the event) in the upper stratosphere of higher latitudes. 
The quite good agreement between modelled and observed 
CIO before the SPEs of January 2005 is expected to im- 
prove modelling results of SPE-induced CIO changes. Be- 
cause CIONO 2 is strictly dependent on CIO, it should be true 
also for CIONO 2 . 

In general, the main objectives of this study are (A) to test 
the ability of WACCM4 to reproduce chlorine changes in- 
duced by January 2005 SPEs; and (B) to validate chemical 
schemes used in current atmospheric models. Other related 
goals are (1) to investigate the temporal evolution of HC1 
changes without the inclusion of the ion chemistry in the 
model; ( 2 ) to check if the favorable conditions of the vor- 
tex, roughly coincident with the terminator, allows a better 
simulation of the changes in CIO and CIONO 2 than in prior 
studies; and (3) to corroborate the hypothesis of chlorine ac- 
tivation under SPEs. 

This paper is divided into eight sections, including the 
introduction. In Sect. 2 we describe the MLS-Aura and 
MIPAS-ENVISAT satellite data used. Section 3 reports an 
overview of the changes in chlorine species induced by the 
SPEs of January 2005 and compares them with the year-to- 
year variability and with the variations caused by the SPEs of 
December 2006. In Sect. 4 we describe the main characteris- 
tics of the version 4 of WACCM adapted to reproduce the im- 
pact of SPEs on the chlorine species. The space weather and 
the meteorological conditions of January 2005 are discussed 
in Sect. 5. In Sect. 6 the spatial distributions and temporal 
development of measured trace gas concentrations are com- 
pared to model calculations and in Sect. 7 these results are 
discussed in the context of previous work. In Sect. 8 we con- 
clude which features in the data are well understood, which 
need further investigation, and how further studies could im- 
prove our current knowledge on SPE induced chlorine chem- 
istry. 

2 Satellite data 

The main features of the satellite instruments and respective 
datasets used in the present study will be discussed in this 
section. 
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2.1 MLS 

The NASA EOS (Earth Observing System) MLS 
(Waters et al., 2006) is one of the four instruments on Aura 
satellite, launched on 15 July 2004 to a sun-synchronous near 
polar orbit with a 13:45 LT (local time) ascending equator- 
crossing time. MLS scans the limb of the Earth in the flight 
forward direction, leading to data coverage from 82° S to 
82° N latitude on every orbit, and measures microwave 
emission in different spectral regions. Vertical profiles are 
measured every 165 km along the suborbital track. Their 
horizontal resolution is 200-300 km along-track and 3-9 km 
across-track. The vertical resolution depends on the species. 
It is typically 3—4 km in troposphere and stratosphere and 
it degrades up to 8 km in the mesosphere. MLS records 
profiles of gas phase chemical abundances, temperature, 
and cloud ice density. Here we will focus mainly on MLS 
Version 3.3 CIO, HC1, HOC1, N 2 O and temperature Level 2 
Data (available at: http://mirador.gsfc.nasa.gov/index.shtml). 
Detailed validation results for the Version 2.2 MLS mea- 
surements are shown in Livesey et al. (2007). Additional 
information on the differences between Version 2.2 and 3.3 
data is reported in Livesey et al. (201 1) and in the individual 
validation papers referenced therein. 

MLS CIO data (Santee et al., 2008) have a precision on 
individual profiles between 0. 1-0.3 ppbv and fairly small 
(lower than 15 %) systematic uncertainties (or good accu- 
racy) in the middle/upper stratosphere. Vertical resolution, 
determined from the full width at half maximum (FWHM) 
of the averaging kernel, is about 3-4.5 km. The useful range 
for a single profile should be between 100 and 1 hPa. MLS 
HC1 data (Lroidevaux et al., 2008) have a precision of 0.2- 
0.7 ppbv in the stratosphere and lower mesosphere and good 
accuracy (10%). The vertical resolution is about 3 km in the 
stratosphere and 4-5 km in the lower mesosphere. The use- 
ful range is between 100 and 0.32 hPa. MLS HOC1 data (see 
Livesey et al., 2011) are quite noisy, with single -profile pre- 
cision of 0. 3-0.4 ppbv, and generally require some averag- 
ing (e.g., weekly zonal mean) to obtain useful precision. In 
general, MLS HOC1 volume mixing ratios (VMRs) seem to 
be consistent with MIPAS HOC1 in the 10 to 2hPa recom- 
mended MLS vertical range. The vertical resolution is about 
6 km. Mixing ratio values where the a priori information has 
a strong influence are flagged with negative precision for 
MLS data, and they should not be used in scientific analyses. 
Although the recommended upper limit for HOC1 profiles is 
2 hPa, HOC1 data used in this study are flagged positively up 
to 1 hPa (comprised); therefore in this region there is not a 
strong influence of the a priori. Moreover, during the space 
weather storm of January 2005, HOC1 VMRs increase by a 
factor of 3 in the upper stratosphere with respect to typical 
pre-SPE conditions. In this way the increased signal-to-noise 
ratio allows us to get information up to 1 hPa. The single- 
profile precision of MLS N 2 O data (Lambert et al., 2007) 
is about 13-25 ppbv (7-38%) and the accuracy is about 3- 
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70ppbv (9-25%) for the pressure range 100-4.6 hPa. The 
vertical resolution determined from LWHM of the averaging 
kernel is 4—6 km. The scientifically useful range of the N 2 O 
data is between 100 to 0.46 hPa. 

Since the magnitude of SPE-induced variability is ex- 
pected to be of the order of (or less than) the individual pro- 
file precision, we will focus on averaged daily data. 

2.2 MIPAS 

MIPAS is a mid-infrared Lourier transform limb emission 
spectrometer designed and operated for measurement of at- 
mospheric trace species from space (Lischer et al., 2008). 
It is part of the instrumentation of the European Environ- 
mental Satellite (ENVISAT) which was launched into its 
sun-synchronous polar orbit of 98.55° N inclination at about 
800 km altitude on 1 March 2002. MIPAS passes the equator 
in a southerly direction at 10:00 a.m. LT 14 to 15 times a day, 
observing the atmosphere during day and night with global 
coverage from pole to pole. The instrument’s field of view is 
30 km in the horizontal direction and approximately 3 km in 
the vertical direction. MIPAS operated during January 2005 
at a reduced spectral resolution of 0.0625 cm -1 (unapodized) 
in terms of full width at half maximum. During this pe- 
riod, MIPAS measurements were performed on 10-13 Jan- 
uary (middle atmosphere observation mode), 16-18 January 
(upper troposphere-lower stratosphere (UTLS) observation 
mode) and 27-28 January (nominal observation mode). Up- 
per atmospheric observation with a scan range of 40-170 km 
in terms of tangent heights, taken on 21-24 and 29-30 Jan- 
uary, are not considered in this study since they provide 
only marginal information on CIONO 2 , the MIPAS key con- 
stituent of this study. The middle atmospheric CIONCL ob- 
servations taken on 10-13 January (i.e., a week before the 
SPE event) turned out to be of limited use in providing un- 
perturbed background CIONO 2 distributions for this study 
because of the high temporal variability of this constituent 
in the vertical range of interest. In the UTLS observation 
mode, MIPAS records a rear- viewing limb sequence of 19 
spectra, corresponding to an along track sampling of approx- 
imately 290 km whereas, in the nominal observation mode, 
it records 27 spectra, corresponding to a sampling of about 
410 km. Tangent heights covered the altitude range from 50 
down to 6 km in the UTLS mode, and 72 to 7 km in the nom- 
inal mode. Trace gas profiles have been retrieved from cal- 
ibrated geolocated limb emission spectra with the scientific 
MIPAS level 2 processor developed and operated by the In- 
stitute of Meteorology and Climate Research (IMK) in Karl- 
sruhe together with the Instituto de Astroflsica de Andalucla 
(IAA) in Granada. The general retrieval strategy, which is a 
constrained multi-parameter non-linear least squares fitting 
of measured and modeled spectra, is described in detail in 
von Clarmann et al. (2003). 

The CIONO 2 data versions used here are 
V4O_CLONO2_100 (UTLS mode) and V4CLCLONO2_201 
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Fig. 1. From top to bottom: time series of daily MLS HC1, HOC1 and CIO at 2hPa (left column) and 1.5 hPa (right column) averaged over 
70-82° N from December to February for 2004-2011 years. The vertical dashed lines show the SPEs occurrence. The estimated precision 
(vertical error bars) for 2004-2005 is reported as a reference. 


(nominal mode), all based on the retrieval setup for the 
original high spectral resolution MIPAS CIONCL mea- 
surements described in Hopfner et al. (2007) (see von 
Clarmann et al. (2009) for the description of the reduced 
spectral resolution CIONO 2 MIPAS retrievals used in this 
study). Differences between the two versions are related 
to the vertical scan range and sampling, implemented in 
the UTLS and nominal observation modes, and do not 
affect noticeably the data characteristics in the altitude 
region of interest (Lopez-Puertas et al., 2005a; Funke et al., 
2011). The single measurement precision ranges from 0.06 
to 0.12ppbv, increasing with altitude. Vertical resolution 
is 5-8 km below 2hPa and 12-14km above. Meaningful 
data is obtained below approximately 0.5 hPa (~52km). 
Unfortunately no MIPAS HOC1 data are available for the 
period under investigation, because since 2005 MIPAS 
measures at reduced spectral resolution not sufficient for 
retrieval of HOC1. 


3 Changes in chlorine species induced by SPEs and 
year-to-year variability 

The descending phase of Solar Cycle 23 was characterized 
by the occurrence of some intense SPEs (e.g., SPEs of Jan- 
uary 2005, May 2005, September 2005 and December 2006). 
In particular, the SPEs of January 2005 and December 2006 
occurred under similar seasonal conditions and impacted 
mainly the atmosphere of the winter hemisphere at Northern 
Polar regions. In order to evaluate the importance of the chlo- 


rine changes induced by the SPEs of January 2005 compared 
with the year-to-year variability and with the effects caused 
by the SPEs of 7 and 13 December 2006, Fig. 1 shows the 
time series of (top to bottom) daily MLS HC1, HOC1 and 
CIO measurements for the winter months (from December 
through February) of the years 2005-201 1 averaged over 70- 
82° N at 2 hPa (left column) and 1 .5 hPa (right column) (note 
that we show two vertical levels fairly close to each other and 
near the upper limit of the MLS data use because chlorine 
changes, caused by the SPEs of December 2006 and by year- 
to-year variability, are visible mainly there). The occurrence 
of the SPEs of January 2005 and December 2006 is shown by 
the vertical black dashed lines. The time series of MLS HC1 
at 2 hPa are illustrated in the upper left panel. Note the strong 
decrease, likely induced by Reaction (R2) and lasting about 
10 days, that occurred after the SPEs of 17 and 20 January 
2005. The HC1 decrease is about 0.3 ppbv, and the minimum 
values occurred the day after the SPE of 20 January. They 
are among the lowest HC1 values ever recorded by MLS at 
2hPa in NH winter during 2004—2011 years. At lower pres- 
sures (up to about 0. 1 hPa; not shown) these features are 
similar, with a significant depletion appearing during SPEs 
and persisting only for a few days. The first proton injection 
of 17 January seems to cause an initial HC1 decrease, over 
which an additional depletion caused by the SPE of 20 Jan- 
uary is imposed (see Sect. 5.1 for the space weather overview 
of January 2005 SPEs). The middle left panel of Fig. 1 shows 
time series of MLS HOC1 measurements as in the upper 
panel. The sudden SPE-induced enhancement of HOC1 due 
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to Reaction (R4) can be easily identified. The HOC1 increase 
(more than 0.2 ppbv) in January 2005 is similar in magnitude 
to the HC1 decrease and suggests SPE-induced chlorine acti- 
vation. These changes are about 3 times higher than the vari- 
ability encountered in the other years. Due to the different 
proton energy spectra of the SPEs of December 2006 com- 
pared with the January 2005 ones, in 2006 the strongest ef- 
fect of the SPE-induced variations is mainly located in the 
upper stratosphere/mesosphere (roughly above 2-3 hPa). In 
contrast, as we will see in the Sect. 6, stratospheric chlo- 
rine changes induced by SPEs of January 2005 are evident 
at least down to 4—5 hPa. Indeed high HOC1 VMRs at 2 hPa 
are also evident during the SPEs of December 2006 (cyan 
line). At this pressure surface, the corresponding HC1 de- 
crease depicted in the upper panel, even though partly hid- 
den by the year-to-year variability, is comparable to SPE- 
induced changes in January 2005. Similar or even more evi- 
dent changes of HC1 and HOC1 for December 2006 and Jan- 
uary 2005 are visible at 1.5 hPa in the upper and middle right 
panels. These important changes, highlighted also in Decem- 
ber 2006, are additional evidence of the presence of chlorine 
activation under SPE occurrence. 

The SPE-induced changes of HC1 were more than 0.2 ppbv 
at 1 .5 hPa, but even larger HC1 decrease took place in Febru- 
ary 2006, 2009 and 2010 corresponding to periods with 
strong winter polar vortices in the upper stratosphere follow- 
ing intense stratospheric sudden warming (SSW) events that 
occurred in January of those years (e.g., Siskind et ah, 2007; 
Manney et al., 2009; Damiani et ah, 2010b). Further HC1 
changes (enhancements) occurred during some SSW events 
(see for example the HC1 peak at 2 hPa in late February 2008 
in the upper-left panel of Fig. 1). Figure 1 shows that HC1 de- 
pletions were coupled with intense HOC1 enhancements only 
during SPEs whereas negligible HOC1 variations occurred 
after or during SSWs when the active chlorine was mainly 
in the form of CIO. On the other hand, no clear evidence of 
CIO changes connected to the SPEs of January 2005 and De- 
cember 2006 is present at 2 hPa. Limited CIO decreases are 
present at 1 .5 hPa even if they cannot be attributable with cer- 
tainty to SPE occurrence. In addition, a sudden spike of CIO 
VMR (about 0.1 ppbv) is visible at 1.5 hPa on 19 January 
2005. We will examine these issues in Sect. 6. 


4 Description of the whole atmosphere community 
climate model (WACCM4) 

In the following, MLS and MIPAS data are compared 
with version 4 of the National Center for Atmospheric 
Research (NCAR) Whole Atmosphere Community Cli- 
mate Model (WACCM4) results (http://www.cesm.ucar.edu/ 
working_groups/WACCM/). An earlier version of the model, 
WACCM3, was used extensively to investigate the atmo- 
spheric effects of SPEs (Jackman et al., 2008, 2009, 2011). 
The present model version, WACCM4, was used in the re- 
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Ion rates (|cm 3 s ’) for Jan. 2005 SPEs 



Fig. 2. Hourly average ion pair production rates (cm -3 s -1 ) as a 
function of time for 16-22 January 2005. The contour levels are 
100, 200, 500. 1000, and 2000cm -3 s -1 . Proton flux data from 
GOES- 1 1 were used to compute ion pair production profiles. 


cent Funke et al. (2011) study of the 2003 Halloween SPEs. 
For this investigation, WACCM4 is forced with meteorolog- 
ical fields from the Modern Era Retrospective-analysis for 
Research and Applications (MERRA, http://gmao.gsfc.nasa. 
gov/research/merra/). The same approach was used in the 
study of the 2003 Halloween SPE (Funke et al., 2011) and 
the study of the 2006 stratospheric sudden warming (Marsh, 
2011). As in those studies, the forcing is achieved by relax- 
ing horizontal winds and temperatures with a time constant 
of approximately 50 h from the surface to 40 km. Above that 
level the forcing is reduced linearly, so that the model is free- 
running between 50 km and the model top. The model has 88 
levels in the vertical, with a resolution in the horizontal of 2.5 
longitude by 1.9 latitude. 

Two sets of simulations were conducted: one with hourly 
energetic proton ionization rates derived from GOES 1 1 flux 
data, and a reference run without ionization. Ionization rates 
are shown in Fig. 2 and discussed in Sect. 5.1. Neither simu- 
lation included highly energetic electron precipitation. This 
is because this study focuses mainly on the middle/upper 
stratosphere where, generally, the electrons-induced ioniza- 
tion is somewhat low. Only for HC1 we have reported data 
up to 0.5 hPa. However, as shown below (see Figs. 4, 5 and 
6), WACCM4 does a good job in reproducing MLS HC1 at 
the upper investigated altitudes, excluding (or at least reduc- 
ing) the possibility of an underestimation of the ionization 
induced by the significant electrons impact on this altitudi- 
nal range. In a future work we will implement the ioniza- 
tion induced by energetic electrons. WACCM4 output clos- 
est to each of the satellite (MLS and MIPAS) geo-locations in 
space and time has been used. It should be noted that in order 
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to obtain more realistic comparisons, MLS or MIPAS aver- 
aging kernels (AK) have been applied to WACCM4 output. 
We did not show the impact of the MLS AK applied to the 
WACCM output because their effect is small. On the other 
hand the effect of the AK is more evident for MIPAS data 
and we have shown this in the right panel of Fig. 7. 


5 Space weather and meteorological conditions during 
January 2005 

Some features of the space weather and the meteorologi- 
cal conditions of the upper stratosphere/lower mesosphere at 
Norther high latitudes during January 2005 are discussed in 
the present section. 

5.1 Space weather 

January 2005 was characterised by low/moderate solar ac- 
tivity between 1 and 13 January and after 23 January, while 
during 15-20 January, the active region 10720 produced five 
powerful solar flares classified as X1.2, X2.6, X3.8, X1.3 
and even X7.1 (see for example The Preliminary Report and 
Forecast of Solar Geophysical Data, called “The Weekly”, 
available at http://www.swpc.noaa.gov/weekly/index.html). 
Three SPEs occurred on 16, 17 and 20 January. In the fol- 
lowing, we will refer exclusively to the SPEs of 17 and 
20 January. The 17 January event was characterised by a 
large injection of solar particles peaking at 5040 pfu (Par- 
ticle Flux Units; 1 pfu= 1 pcm -2 sr -1 s -1 ) at the lOMeV 
channel of the National Oceanic and Atmospheric Admin- 
istration (NOAA) Geostationary Operational Environmental 
Satellite (GOES) while the 10 MeV peak of the 20 January 
event was 1860 pfu. Nevertheless, 20 January exhibited the 
hardest and most energetic particle event of the solar activ- 
ity cycle 23 (see for instance Damiani et al., 2008). Indeed, 
on 20 January the peak for the 100 MeV channel of GOES 
was 652 pfu (we recall that during the Halloween Storm the 
peak flux at the 100 MeV channel was 155 pfu; a list of SPEs 
affecting the Earth’s environment is provided by NOAA at 
http ://w ww. s wpc . noaa.gov/ftpdir/indices/SPE.txt) . 

The different energy spectra of the two events lead to im- 
pacts at quite different layers of the polar atmosphere. In par- 
ticular, the SPE of 20 January is expected to influence the 
atmosphere to a deeper level than the 17 January event due 
to its larger flux of protons at high energy. 

It is important to state that the ion production rate was 
mostly smaller in both 17 and 20 January 2005 compared 
with the SPEs of October 2003 and July 2000. For example, 
in the mesosphere, roughly between 1 and 0.01 hPa (about 
50-80 km), the daily average peak of ion pair production 
rates was more than 5000 cm -3 s - 1 on 29 October 2003 and 
on 15 July 2000, whereas it was about 1000 cm -3 s -1 on 
17 January 2005 (Jackman et al., 2008, 201 1). The largest ion 
production rate in 2005 was located at higher altitudes com- 
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pared with 2003. That influenced the mesospheric NO x pro- 
duction by moving the SPE-induced peak of NO x in January 
2005 toward higher altitudes compared with 2003. Moreover, 
since the SPE of 20 January was characterised by a very short 
proton injection, it is important to take into account a tempo- 
ral resolution of the ion production rate finer than that used in 
some prior studies (e.g., Jackman et al., 2008, 2011). There- 
fore, Fig. 2 shows the hourly average ionization rate, from 
100 hPa to 0.001 hPa, as a function of time for 16-22 Jan- 
uary 2005. The proton flux data from GOES-11 were used 
to compute ion pair production profiles used in the model 
(see Jackman et al., 1980, for the methodology). We can 
note some structure that is not resolvable when looking at the 
daily average ionization rate (compare this figure with Fig. 1 
of Jackman et al., 2011). The most intense ionization, with 
peak greater than 2000 cm -3 s -1 , is above 1 hPa on 17 Jan- 
uary whereas it is below 1 hPa on 20 January. However high 
ionization rates last only a few hours on 20 January. 

5.2 Meteorological conditions 

The absence of strong temperature increases and associated 
reversals of the zonal winds indicates that the NH winter of 
2004/05 was without major sudden stratospheric warmings 
(e.g., Chshyolkova et al., 2007; Xu et al., 2009; Damiani 
et al., 2010b). However, the winter of 2004/05 was partic- 
ularly cold and dynamically active, especially in the lower 
stratosphere (Manney et al., 2006; WMO, 2007; El Amraoui 
et al., 2008). The temperatures exhibit rapid fluctuations in 
some periods, for example at the beginning and at the end 
of January. While the polar vortex reached its strongest state 
and was centred on the pole in December, it was elongated, 
slightly shifted off the pole and had a decreased extension in 
the lower stratosphere during these disturbed periods of Jan- 
uary (Chshyolkova et al., 2007). Because one of these periods 
occurred a few days after the SPEs occurrence, this anoma- 
lous dynamical vortex behavior likely weakened the chemi- 
cal impact of the SPEs. 

Since meridional transport and mixing can play a key role 
in the magnitude and spatial distribution of SPE-induced 
variations. Fig. 3 reports the zonal mean geographic dis- 
tribution of N 2 O in four periods (from top to bottom: 12- 
lb January, 17-21 January, 22-26 January and 27-30 Jan- 
uary) during the second half of January 2005 as recorded by 
MLS (left panels) and simulated by WACCM4 (right pan- 
els). Due to the strong variation of N 2 O with altitude, we 
show the N 2 O relative anomaly, with respect to the merid- 
ional average (40-82° N), expressed in percent. Lower abun- 
dances of N 2 O, surrounded by a pronounced latitudinal gra- 
dient, indicate the position of the polar vortex. During the 
period of the SPEs (17-21 January) MLS shows that the 
vortex is well defined in the stratosphere up to 2-3 hPa, 
whereas above this layer, a smaller latitudinal gradient is 
indicative of more intense latitudinal redistribution. This is 
roughly in agreement with the analysis of the NH polar 
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Fig. 3. Zonal mean relative anomaly of N 2 O, with respect to the average (40-82° N), expressed in percent. From top to bottom: zonal mean 
anomaly of N 2 O before (12-16 January), during (17-21 January) and after (22-26 January and 27-30 January) the SPEs of 17 and 20 January 
2005. Left panels: MLS data. Right panels: WACCM4 result. 


vortex dynamics of Verronen et al. (2011a) showing increas- 
ing horizontal transport and mixing over the altitude range 
from the lower (~ 1 hPa) to middle mesosphere during the 
January 2005 SPEs (however note that the small latitudi- 
nal gradient of N 2 O could be indicative also of a well de- 
fined but displaced vortex). The area 70-82° N, to which we 


will refer in the paper later on (see, for example. Fig. 4), 
seems to be well inside the vortex (with the consequent 
weak redistribution) only up to about 3 hPa. Above this pres- 
sure level, the vortex gets weaker particularly during the 
following days (lower panels). Therefore SPE-induced ef- 
fects should be longer lasting in the middle stratosphere 
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compared with the effects in the upper stratosphere/lower 
mesosphere. WACCM4 reproduces the MLS N 2 O pattern 
quite well but simulates weaker latitudinal redistribution in 
the mesosphere. The temporal evolutions of daily mean 
changes of temperature and N 2 O with respect to January av- 
erage, averaged over 70-82° N and 58-70° N, are reported in 
the upper panels of Figs. 4 and 5. The patterns of the MLS 
temperature and N 2 O changes (left panels) are quite similar 
to each other and show two major variations at the begin- 
ning and the end of January connected to the disturbed po- 
lar vortex periods. In general higher temperatures and N 2 O 
values are associated with entrainment of external vortex 
air into the 70-82° N region. Thus they can indicate intru- 
sions or vortex displacement. However, the isolines of the 
absolute N 2 O abundances (not shown) suggest continuing 
confined air descent in the middle stratosphere throughout 
January. Therefore in this region the vortex mixing barrier 
should remain quite strong at this time according to Chshy- 
olkova et al. (2007). More frequent variations of temperature 
and N 2 O occur in the upper stratosphere, so, as already dis- 
cussed, this region of the vortex is less stable compared with 
lower altitudes. Further, in Fig. 4 we observe an increase of 
N 2 O values in the upper stratosphere starting around 18 Jan- 
uary and becoming more pronounced towards the end of the 
month. Inspection of the MLS N 2 O data maps (not shown) 
suggests that, although the polar vortex remains substantially 
well defined in January, it is elongated and shifted away from 
the North Pole on many days. In particular, after 27 Jan- 
uary, a considerable part of the region inside 70° N is oc- 
cupied by extra-vortex air and obviously that influences the 
daily means. The WACCM4 simulations follow the changes 
in MLS temperature and N 2 O quite well, although some dis- 
crepancies in their magnitude are visible in the second half 
of January. Finally, zonal mean changes of temperature and 
N 2 O on 17-21 January with respect to 12-16 January at 40- 
90° N are reported in the upper panels of Fig. 6. Temperature 
and N 2 O patterns are consistent with the general behaviour 
of the vortex already depicted, hence we observe its intensi- 
fication in the middle stratosphere and its weakening in the 
upper stratosphere. 

6 Changes in chlorine species: satellite data vs. 

WACCM4 

To provide the best representation not only for SPEs but also 
for the meteorological evolution. Fig. 4 shows the temporal 
evolution of daily mean changes of HC1, HOC1, CIO with re- 
spect to January average, averaged over 70-82° N, therefore 
under almost complete nighttime conditions (see Sect. 5.2 
for the description of temperature and N 2 O panels). The loss 
in MLS HC1 (left panel) occurred after 17 January and af- 
fects nearly all investigated pressures. The strongest decrease 
(about 0.25 ppbv) occurs uniformly in a rather wide region 
ranging from 1.5 to 6hPa on 21-22 January. This decrease 
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lasts until the end of January slightly descending to around 
6-8 hPa. WACCM4 simulations of HC1 (central panel) show 
a quite good agreement in the lower mesosphere but not in 
the stratosphere. The HC1 decrease below 2hPa, well evi- 
dent in MLS data, is not reproduced by WACCM4 results. 
The right panels of Fig. 4 (WACCM p-b) show the differ- 
ence between the perturbed (with SPEs) and the base (with- 
out SPEs) WACCM4 simulations. Comparing WACCM4 re- 
sults with WACCM(p-b) shows that, in WACCM4 results, 
the HC1 decrease caused by the SPEs between 2 and 4-5 hPa 
is masked by background variability, and the HC1 depletion 
below 4-5 hPa could be only slightly influenced by the SPE 
impact. The MLS HOC1 temporal evolution shows an in- 
crease of more than 0.2 ppbv on 18-21 January 2005. After 
the first enhancement caused by the SPE of 17 January (more 
than 0.15 ppbv), a further HOC1 increase is induced by the 
SPE of 20 January (see Fig. 4) and HOC1 disturbances are 
slightly moved downward because of the harder spectrum of 
the second SPE. This enhancement is clearly evident above 
4-5 hPa and lasts up to about 25 January. Both the pattern 
and the magnitude of these changes are well reproduced by 
WACCM4. 

The lower left panel of Fig. 4 shows the temporal evolution 
of the MLS CIO. CIO presents a weak increasing trend in the 
region 3-5 hPa, that reaches the highest values after 20 Jan- 
uary, and a decrease in the upper stratosphere (1-2 hPa) 
which is strongest around 18 January. These changes seem 
to be mainly related to dynamical variability. The decrease 
of CIO starts before the occurrence of the SPEs and there- 
fore it masks the influence of the solar forcing. However, a 
further depletion in CIO occurs on 17 and 18 January. Then, 
on 19 January a slight CIO rise is quite evident, followed 
by another depletion on 20 and 21 January. This behaviour is 
compatible with a reduction of CIO via Reaction (R4) and the 
consequent HOC1 formation. The WACCM4 result (central 
panels) reproduces the general features of MLS CIO. How- 
ever, the onset of the particle-induced CIO depletion is more 
clearly visible since the observed decreasing trend before the 
SPEs is less pronounced in the simulation. The rather con- 
stant values of the WACCM(p-b) differences (right panel) 
below ~4hPa suggest that the observed and modelled CIO 
increase is not connected to the SPE. 

Figure 5 shows the temporal evolutions as in Fig. 4 but 
with values averaged over 58-70° N, and therefore roughly 
under twilight conditions. In this geographic region MLS 
HC1 still undergoes a significant decrease (up to about 
0.2 ppbv) that is mostly evident above 1 hPa and around 
5 hPa. As for the 70-82° N region (Fig. 4), WACCM4 does 
not reproduce the HC1 depletion below about 2 hPa. More- 
over, unlike the Fig. 4, it clearly underestimates the meso- 
spheric HC1 depletion. From the comparison of WACCM4 
results and WACCM(p-b) differences it follows that the sim- 
ulated HC1 decrease below 4 hPa is likely not due to SPE. 
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Fig. 4. From top to bottom: temporal development of daily MLS T, N 2 O, HC1, HOCI, CIO during January 2005, averaged over 70-82° N, 
relative to background conditions (January average). Left panels: MLS data. Middle panels: WACCM4 result. Right panels: the model SPE 
response (perturbed - base). 
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It is interesting to examine the short-term variations of 
HOC1 and CIO as observed by MLS during the SPEs oc- 
currence. At 58-70° N, under greater solar illumination with 
respect to the 70-82° N region, an HOC1 enhancement, with 
peaks of about 0.1 ppbv on 18, 20 and 22 January, is evident 
around 2-3 hPa. CIO changes are also evident. In general, a 
CIO decrease during SPE-influenced days, as seen in Fig. 4, 
is not present. On the contrary, a blurred CIO enhancement 
takes place around 1 hPa. Moreover, CIO increases are quite 
evident on 19 and 21 January, i.e. not in conjunction with 
but the day after the HOC1 peaks (note that the CIO peak 
of 19 January is evident also at nighttime, see Figs. 1 and 4). 
This suggests that the SPE-generated HOC1 is photolyzed af- 
ter the SPEs, leading to CIO formation via Reaction (R3). 
In addition, somewhat smaller CIO increases seem to be 
present also around 4-6 hPa, roughly coincident with the re- 
gion where stratospheric HC1 decrease occurred. WACCM4 
reproduces very well the HOC1 increases whereas it does not 
show any CIO enhancement. On the contrary it shows a very 
small CIO decrease (as at nighttime), coupled with the SPEs 
occurrence, that is not present in the MLS data. 

Figure 6 illustrates the zonal mean changes of (from top 
to bottom) HC1, HOC1, CIO on 17-21 January 2005 with 
respect to 12-16 January 2005 at 40-90° N (see Sect. 5.2 
for the description of temperature and N 2 O panels). The ge- 
ographic distribution of MLS HC1 and HOC1 changes (left 
panels) shows that the perturbations above 2 hPa extend well 
also into mid-latitudes. This feature is less pronounced in 
WACCM4 (central panels), most probably related to the un- 
derestimation of the latitudinal redistribution at higher alti- 
tudes. Also below 4hPa negative changes of MLS HC1 ex- 
tend towards lower latitudes, however they are not connected 
to SPEs occurrence. N 2 O changes at 50-70° N (evident in 
the above panels) are indicative of the different behavior of 
the polar vortex in these two periods. It was roughly circu- 
lar in shape and centered over the pole on 17-21 January 
whereas it was elongated and slightly shifted off the pole on 
the days before. Therefore, taking into account differences 
in HC1 abundance between inside and outside of the polar 
vortex, dynamics is likely the main reason of the negative 
changes in HC1 at 50-70° N visible in both observations and 
modelling results. 

Both HC1 and HOC1 show strongest changes at the high- 
est latitudes, where the photolysis is lowest and the vortex 
is more isolated. Because of the limited ionization at lower 
altitudes, only HC1 changes above about 4-5 hPa are likely 
due to SPEs occurrence (see also the WACCM(p-b) in the 
right panel of Fig. 6). On the other hand, the evident model- 
data HC1 discrepancies below 4—5 hPa inside the vortex could 
hardly be attributed to unrealistic model dynamics, given the 
good agreement of observed and modeled tracer (and tem- 
perature) evolutions, there. So we cannot exclude that some 
small SPE-induced HC1 depletion could occur also in this re- 
gion. However, it could be possible that this inconsistency is 
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caused by a different chemical response to dynamical back- 
ground variations in the model. 

The zonal mean change of CIO (lower panels of Fig. 6) 
shows a well defined pattern in both MLS data and 
WACCM4 results. The observed CIO depletion of about 
0. 1 ppbv is present only at high latitudes whereas a CIO en- 
hancement of a similar magnitude is visible at about 60- 
70° N. WACCM4 results underestimate the response at lower 
latitudes and overestimate the vertical extent of the CIO de- 
crease at high latitudes. The different response of CIO to 
SPEs as a function of latitude seems to be related to the pres- 
ence/absence of solar illumination. Therefore, HOC1 is the 
main active chlorine species during nighttime, whereas both 
CIO and HOC1 show perturbations, probably due to SPEs, 
in twilight regions. However, WACCM(p-b) differences sug- 
gest that the CIO increase at lower latitudes is not linked to 
SPE occurrence (see also the WACCM(p-b) of Fig. 5). We 
will reexamine these features in the discussion section. 

MIPAS and WACCM4 CIONO 2 zonal mean changes with 
respect to 16-17 January 2005 are shown in Fig. 7 for 18 and 
27 January 2005. In the upper panels of Fig. 7 we can see 
that the zonal mean change on the day after the first SPE does 
not indicate a SPE-induced CIONO 2 response; observed and 
simulated small changes on the order of 0.03-0.04 ppbv are 
background variability (see WACCM(p-b) differences). The 
lower panels show that the zonal mean change on 27 Jan- 
uary is characterised by large CIONCL changes (~0.2ppbv) 
around 1 hPa and by a smaller variation (~0.07 ppbv) at 
about 3 hPa. WACCM4 results reproduce well the observa- 
tions in terms of magnitude and altitude of these peaks. In 
addition, WACCM(p-b) differences suggest that the CIONO 2 
changes around 3 hPa are mainly attributable to SPE occur- 
rence, whereas the enhancement at 1 hPa is induced by both 
background variability and SPE occurrence, with the latter 
contributing ~0.1 ppbv mainly at latitudes higher than 70° N 
and between 1 and 2 hPa. 


7 Discussion 

The analysis presented above highlighted the presence of 
chlorine activation under SPE conditions although some as- 
pects still remain unclear. Therefore, it is worthwhile to dis- 
cuss these results in the light of past findings. Some impor- 
tant features distinguish the January 2005 SPEs from other 
events (e.g., SPEs of July 2000 and October-November 2003) 
investigated in the context of their impact on chlorine. First, 
the polar night region was much more extended in January 
2005 compared with October 2003 and July 2000. The well 
established vortex coincided with the terminator in 2005, 
while in 2003, the early winter vortex was partially illumi- 
nated, and the SPE of July 2000 was investigated only in the 
summer hemisphere. It is also important to emphasise the 
smaller ionization rates modeled for January 2005 compared 
with other events (see Sect. 5.1), which caused a smaller 
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Fig. 6. From top to bottom: zonal mean change (17-21 January 2005 with respect to 12-16 January 2005) of MLS T, N 2 O, HC1, HOC1, CIO 
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Fig. 7. CIONO 2 zonal mean change on 18 (upper panels) and 27 January 2005 (lower panels) with respect to 16-17 January 2005 at 40- 
90° N. From left to right: MIPAS data, WACCM4 result with MIPAS averaging kernels applied, the model SPE response (perturbed - base), 
WACCM4 result with MIPAS averaging kernels applied - the original model result. 


production of NO x at upper altitudes in 2005 compared to 
2003. As we will see in the following, these elements are 
important in order to explain the impact of the investigated 
events on the atmospheric composition. 

7.1 HC1 changes 

To our knowledge, only Winkler et al. (2009) have carefully 
examined the variation of HC1 under SPE conditions. How- 
ever, the different time of the year of the examined SPE 
and the analysis performed by means of solar occultation 
data and a 2-D chemistry model (with ion chemistry) pre- 
vent a full comparison of their results with ours. Some dif- 
ferences in the vertical distribution of the HC1 depletion, 
caused by the different SPEs, are evident. HALOE data in 
2000, recorded in the illuminated NH, show variability be- 
tween 65 km (~0. 1 hPa) and 35 km (~5 hPa), with the max- 
imum stratospheric HC1 depletion (about 0.4 ppbv) at 50 km 
and a roughly constant reduction of this enhancement down 
to 35 km, where the HC1 changes are small (~0.05 ppbv). In 
contrast, MLS data in January 2005 show the maximum HC1 
loss (of less than 0.25 ppbv) uniformly distributed in a rather 
wide region ranging from 1 .5 to 6 hPa. These differences are 
mainly related to the stratospheric background variations in 
January 2005, to the different atmospheric seasonal condi- 


tions at Northern Polar regions in July and January, and to 
the occurrence of two successive SPEs in January 2005. 

In fact, it is important to identify the distinct contributions 
arising from different sources to the HC1 decrease in January 
2005. In January 2005 we have the overlap of the effects of 
two SPEs, with the second one, characterized by a harder 
spectrum, able to impact the already low stratospheric HC1 
values persisting longer at nighttime (Fig. 4). On the con- 
trary, only one SPE (even if more intense) occurred in July 
2000. The diverse proton energy spectra during the two SPEs 
in January 2005 is reflected in the shape of the HC1 decrease. 
The first depletion (~0.2ppbv), due to 17 January SPE, af- 
fected mesospheric altitudes on 18 January. In this case, the 
location of the maximum HC1 depletion (around 1 hPa) is 
roughly in agreement with previous results. Then, a second 
and more intense stratospheric decrease, superimposed over 
the previous one, occurred after the 20 January SPE. This 
strongest decrease occurs rather uniformly in the upper and 
middle stratosphere on 21-22 January and it lasts until the 
end of January around 6-8 hPa. However, because of the lack 
of significant ionization at these altitudes, this long lasting 
tail of HC1 decrease seems not to be connected to the SPE oc- 
currence. The expected low SPE-induced impact below about 
4 hPa is likely superimposed over background variations (see 
Sect. 6). It results in a stronger than expected HC1 decrease 
with a maximum around 5 hPa on 21-22 January and in the 
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evident HC1 tail persisting until the end of January. Hence, 
these model-data discrepancies, mainly present at about 5- 
8 hPa (see Figs. 4 and 5), are mostly due to background at- 
mospheric variability. Their investigation is beyond the scope 
of this study and they will require a specific analysis. 

An additional element that could play a role in explaining 
these differences in the HC1 decrease is related to the differ- 
ent background of OH VMRs in January and July in Northern 
Polar regions. The OH VMRs in the stratosphere and meso- 
sphere depend mainly on the Solar Zenith Angle (SZA). In 
the upper stratosphere, between 7 hPa and 1 hPa, the typi- 
cal daytime OH abundance is about 0. 1-0.6 ppbv, with lower 
values at lower altitudes, whereas during the night the OH 
VMRs are uniformly very close to zero. In order to have ap- 
preciable HC1 changes, the SPE-induced OH increase has to 
be, at least, of the same order of magnitude (or larger) as the 
OH background. Therefore differences in the vertical distri- 
bution of the HC1 depletion, caused by the SPEs of July 2000 
and January 2005, had likely been influenced by differences 
in the background of OH VMRs. 

A further difference between the SPEs of January 2005 
and July 2000 is the duration of the mesospheric HC1 de- 
crease. Indeed, HC1 recovered in a couple of days in July 
2000, whereas the depletion lasted more than one week in 
January 2005. Obviously the occurrence of two successive 
SPEs in January 2005 is the first reason of the correspon- 
dent longer HC1 depletion. Moreover, the confinement main- 
tained by the polar vortex in January 2005 and the absence 
of photolysis at high latitudes, that likely prevented a quick 
shift towards HC1 from the SPE-induced active forms (e.g., 
the destruction of HOC1 by photolysis), are additional causes 
compared to the polar summer conditions in July 2000. 

WACCM4 reproduces quite well the HC1 depletion in the 
lower mesosphere (about 0.2 ppbv) at nighttime on 18 Jan- 
uary. Also at the upper mesospheric altitudes where the HC1 
depletion is even more intense (not shown), WACCM4 re- 
sults are in fair agreement with experimental data without the 
necessity of including ion chemistry. However, WACCM4 
underestimates the mesospheric HC1 response at 58-70° N 
(see Fig. 5). This could be related to the underestimation of 
the redistribution of the bulk of depleted HC1 molecules. On 
the other hand WACCM4 slightly overestimates the HC1 de- 
pletion in the nighttime mesosphere at 70-82° N (see Fig. 4) 
on 21 January. Since Verronen et al. (2011a) showed a strong 
mesospheric horizontal transport on 21 January, these model- 
data differences are a further corroboration of the imper- 
fect reproduction of the mesospheric horizontal mixing by 
WACCM4 during this period. 

In the upper stratosphere (roughly above 2-3 hPa), where 
HC1 changes are related to the SPEs, the model tends to un- 
derestimate the observed depletion (see Figs. 4 and 5). The 
lack of HC1 loss (or its underestimation) in the WACCM4 re- 
sults in this region is in line with the difficulties encountered 
by Winkler et al. (2011) when modeling the stratospheric 
HC1 response during the July 2000 event. The cause is not 
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easily understandable. Possible reasons that could contribute 
to explain such difference are: 

1. SPE-induced ionization included in the model is too 
small, either due to an underestimation of proton fluxes 
or due to the neglection of electron-induced ioniza- 
tion (Verronen et al., 2011b). However, MLS OH is re- 
produced by WACCM4 quite well down to 3 hPa (not 
shown), so this possibility appears to be unlikely. 

2. The ion chemistry in the model is neglected. 

3. The uncertainty of the rate constant of Reaction (R5) 
included in the model with the JPL rate constant (Sander 
et al., 2006). It is connected to the deactivation of CIO 
into the reservoir HC1 but probably could play only a 
limited role in this discrepancy. 

This model-data discrepancy suggests that deficiencies in 
chemical schemes involving chlorine species in the upper 
stratosphere could be present. Further analyses will be nec- 
essary in order to explain the model-data discrepancies. 

7.2 CIO and HOC1 changes 

It is also interesting to examine in detail the behaviour of 
CIO and HOC1. von Clarmann et al. (2005) showed an in- 
crease of MIPAS CIO VMRs immediately after the SPE of 
September-October 2003. A recent analysis based on repro- 
cessed IMK/IAA MIPAS data (Funke et al., 2011) clarified 
that MIPAS recorded increased CIO (about 0.1 ppbv) mainly 
outside the polar night region, whereas a CIO decrease (about 
0.2 ppbv) occurred in the polar night region. This is consis- 
tent with the MLS results presented in Fig. 6 where the CIO 
decrease is roughly confined to the polar night sector and in- 
creased CIO is evident at the terminator (in particular above 
2 hPa at 60-70° N). However, WACCM4 simulations indi- 
cate that only the CIO decrease is caused by SPEs. On the 
other hand, the strong HOC1 enhancement is well reproduced 
by WACCM4 at nighttime conditions. 

During the daytime, HOC1 is quickly destroyed by photol- 
ysis 

HOC1 + h v (X < 503 nm) -> Cl + OH (R6) 

and by reaction with OH (or atomic oxygen): 

OH + HOCl-> H 2 0 + C10. (R7) 

Funke et al. (2011) suggested that, because OH is the dom- 
inant HO x constituent during daytime (at least in the upper 
stratosphere; below roughly 5 hPa HO 2 starts to be compara- 
ble to OH and gets even larger than OH at lower altitudes) 
and Reactions (R6) and (R7) are the major daytime pro- 
cesses, active chlorine is mainly in the form of CIO under 
these conditions. The CIO enhancements in MIPAS data out- 
side the polar region during the Halloween storm could then 
be related to SPE-induced chlorine activation if the produc- 
tion of HO x was above the typical HO x background levels. 
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In January 2005 MLS data indicate that the SPE-induced OH 
rise is comparable with background levels (about 0.05 ppbv) 
at 1 hPa and 65-70° N, coinciding with the altitude of the 
CIO peak (see left panel of Fig. 6), whereas there are not 
appreciable changes at lower altitudes. 

Figure 6 shows that the region of low HC1 values ex- 
tends towards low latitudes in the upper (~1 hPa) and mid- 
dle stratosphere (~7 hPa), whereas it remains more confined 
inside high latitudes around 3 hPa. The MLS CIO structure 
in the same region (roughly 60-70° N) presents higher CIO 
values at the upper altitudes, where there is availability of 
active chlorine molecules arising from the HC1 depletion via 
SPE-induced OH, whereas little change was recorded around 
3 hPa. Further, HOC1 enhancements are present also at 60- 
70° N above about 3 hPa (Fig. 6). Increased production of 
HOC1 in the middle stratosphere, induced by accelerated odd 
hydrogen chemistry (via Reaction R4), has been reported by 
von Clarmann et al. (2009) in association with a displace- 
ment of the polar vortex towards lower latitudes during a 
SSW event. SPE occurrence represents a further example 
of accelerated odd hydrogen chemistry as indicated by both 
HOC1 and CIO increases between 58-70° N (see Fig. 5). 
Comparing the temporal evolution of CIO and HOC1 zonal 
means, it follows that CIO peaks seem to temporally follow 
HOC1 changes. This suggests that SPE-induced HOC1 in- 
creases may have been photolyzed leading to the shift from 
HOC1 to CIO. Therefore, clues for a deactivation later af- 
ter the SPEs of January 2005 seem to be present. However, 
WACCM(p-b) of Figs. 5 and 6 does not show any SPE- 
induced CIO enhancements whereas it reproduces well the 
HOC1 changes also outside the polar night terminator. The 
dynamically active polar vortex could also have contributed 
to complicate the situation by masking and mixing chlorine 
changes. 

Considerable uncertainty persists in laboratory measure- 
ments of the rate constant for the Reaction (R4). Some stud- 
ies (e.g., Kovalenko et al., 2007; von Clarmann et al., 2009) 
pointed out the necessity of increasing it towards the value 
from Stimpfle et al. (1979) that is about a factor-of-two faster 
than the rate constant (Sander et al., 2006) used in WACCM4. 
Recently, the new JPL recommendation (Sander et al., 201 1) 
was updated towards a higher value. However this remains 
an open issue and the possibility of increasing the agreement 
between observed data and model results by uploading the 
rate constant should be taken into account. 

The SPE-induced HOC1 changes recorded by MLS in Jan- 
uary 2005 are similar to the HOC1 variations recorded by MI- 
PAS during the Halloween storm (about 0.25 ppbv). Even if 
we consider only the enhancement induced by the first SPE 
on 18 January (close to 0.2 ppbv at 1.5 hPa), therefore ignor- 
ing the additional increase due to the second SPE, the de- 
pletion still remains comparable. Because of the lower ion- 
ization of January 2005 (about 4-5 times less, see Sect. 5.1) 
this seems difficult to explain. The amount of SPE-induced 
HOC1 enhancement is linked to the abundance of HCL that 
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depends on the ionization; however also the seasonal condi- 
tions could play an important role. Since, under nighttime 
conditions, HO x is in steady state also during a SPE, its 
abundance is proportional to ionization. The HO x produc- 
tion induced by SPEs occurs mainly in the form of OH and 
H (Solomon et al., 1981). Then, OH can react with ozone 
and H, in presence of a third body, with oxygen molecules, 
both leading to HO 2 formation. As the production of HO 2 
depends on the density, it follows that the HO 2 /OH ratio in- 
creases from the mesosphere down to the stratosphere under 
SPE conditions. Moreover, since the background of OH is 
higher than HO 2 in the daytime upper stratosphere and meso- 
sphere (i.e., the region affected by SPE), compared with its 
very low abundances at nighttime, this ratio is more elevated 
at nighttime. On the other hand, HO x losses are dominated 
by reactions between OH and HO 2 . In this way, the HO 2 /OH 
ratio, i.e. HO x partitioning, influences the HOC1 formation 
since it controls the availability of HO 2 . Therefore, the higher 
SPE-induced HO 2 /OH ratio in January 2005, due to the al- 
most complete nighttime conditions, contributes to an HOC1 
enhancement similar to that of 2003 despite the lower ioniza- 
tion and total HO x production. 

The recent studies of Jackman et al. (2011) and Funke 
et al. (2011) highlighted that models tend to overestimate the 
SPE-induced production of HO 2 and H 2 O 2 and, hence, the 
HO 2 /OH ratio under SPE conditions. As discussed, this ratio 
plays an important role in the partitioning of chlorine species. 
The MLS HO x data could help us to constrain the HO 2 /OH 
ratio. However the presence of a bias in MLS HO 2 (Pickett 
et al., 2008), which normally requires taking day-night dif- 
ferences to be removed, prevents an accurate analysis (SPE- 
induced HO 2 variations are evident mainly under nighttime 
conditions). The biases are typically small compared with 
SPE effects. Using differences of the nighttime data between 
two close time periods (to minimize the bias effect) it is pos- 
sible to highlight SPE-induced HO 2 changes (Jackman et al., 
2011) although in this way the exact HO 2 /OH ratio can be 
affected. 

7.3 CIONO 2 changes 

Also the behaviour of CIONO 2 deserves special attention. 
Figure 7 shows that CIONCL changes are not present on 
the day after the first SPE. On 27 January we note an en- 
hancement of ~0.2 ppbv at about 1 hPa. The work of Funke 
et al. (2011), who use a reprocessed version of the CIONO 2 
data set presented by von Clarmann et al. (2005), shows the 
SPE-induced CIONO 2 peak (~0.4ppbv) at about 3hPa in 
October-November 2003. Either the altitude or the magni- 
tude is somewhat different compared with the peak value ob- 
served on 27 January 2005. The CIONO 2 formation depends 
on the availability of NO 2 and CIO, and on the density. In 
January 2005 the SPE-induced production of odd nitrogen 
was considerably lower compared to the October-November 
2003 events. Jackman et al. (2008) estimated that the NO y 
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production in 2005 was ~4 times smaller than in 2003, which 
could lead to a lower production of ClONCh in January 2005. 
Moreover, we could expect a CIONO 2 enhancement at some- 
what lower altitudes compared with the altitudes of the SPE- 
induced NO 2 layer (because of the pressure dependence of 
CIONO 2 formation). Because of the different distribution of 
the ionization (see Sect. 5.1), the SPE-induced NO 2 forma- 
tion in January 2005 should occur mainly at higher altitudes 
compared with 2003. Hence, it is likely that the CIONO 2 in- 
crease (~0.2ppbv) of 27 January at about 1 hPa (Fig. 6 ) is 
due, at least partly, to the SPEs. This seems compatible with 
the incomplete recovery of HC1 around 1 hPa lasting until 
the end of January (see Figs. 4 and 5) caused by an initial 
chlorine activation (HC1 -» CIO) and a further possible de- 
activation into the CIONO 2 reservoir. The elevated altitude of 
HNO 3 enhancements that occurred during the SPEs of Jan- 
uary 2005 and recently reported by Verronen et al. (2011a) 
(with peak around 60 km, so at altitudes higher than previ- 
ous observed HNO 3 variations caused by SPEs (e.g., Lopez- 
Puertas et al., 2005a) is a further element corroborating our 
finding of SPE-induced CIONO 2 changes at elevated alti- 
tudes. 

It is puzzling that CIO and CIONO 2 are quite well repro- 
duced by WACCM4 in January 2005 but not in October- 
November 2003. Funke et al. (201 1) argued that in 2003 the 
background CIO was underestimated by models, leading to 
modelled CIO changes smaller than observed. In contrast, in 
2005 WACCM4 reproduces quite well the background CIO 
VMRs, and hence also the modelled SPE-induced changes 
are closer to observed data. A possible explanation could be 
that the chemical schemes of CIO formation/loss under twi- 
light conditions are not well represented in the models. An- 
other possible explanation could be that models underesti- 
mate the actual isentropic mixing within the vortex. In 2003, 
a weaker mixing of “illuminated” and “dark” air masses in- 
side the partially illuminated vortex could lead to lower mod- 
elled CIO abundances in the polar regions. In contrast, with 
a well established vortex coincident with the terminator as 
in 2005 the modelled results should be closer to the satellite 
data. The lower CIO abundance in the models during Octo- 
ber 2003 is the main reason for the model underestimation 
of CIONO 2 . In contrast, in January 2005 the modelled CIO 
VMRs are similar to MLS CIO values, hence also WACCM4 
CIONO 2 output reproduces the MIPAS data well. 


8 Conclusions 

For the first time the main contributors to the Cl y family 
recorded by two different satellites have been shown and 
compared with a state-of-the art chemistry climate model un- 
der SPE conditions. This has been done taking into account 
also HC1, the most important chlorine reservoir, to which 
many uncertainties (confirmed by the present study), linked 
to its response to SPEs, are related. MLS measurements of 
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chlorine species (i.e., HC1, CIO and HOC1), N 2 O, and tem- 
perature together with MIPAS CIONO 2 recorded during and 
after the SPEs of 17 and 20 January 2005 have been analysed 
and modelled with version 4 of the Whole Atmosphere Com- 
munity Climate Model (WACCM4). SPE-induced chlorine 
repartitioning (HC1 — > CIO -> HOC1 and, after the events, 
CIONO 2 ) has been identified. In the observations, HC1 de- 
crease is present at almost all the investigated altitudes (i.e., 
10-0.5 hPa). The strongest decrease (about 0.25 ppbv) occurs 
uniformly in a rather wide region ranging from 1.5 to 6 hPa 
on 21-22 January and it recovers in a week in the upper 
stratosphere. However below 4 hPa, this decrease is unlikely 
to be related to the SPE. Instead, it seems to be related to 
background variations. 

CIO abundances did not undergo very large changes. How- 
ever, a CIO decrease of up to 0. 1 ppbv has been identified 
in the upper stratosphere at high latitudes, whereas a CIO 
increase of similar magnitude has been found close to the 
terminator. An HOC1 enhancement of more than 0.25 ppbv 
occurred roughly in the same region where the CIO deple- 
tion took place. This confirms that HOC1 is the main active 
chlorine species produced by proton forcing under nighttime 
conditions. Lower HOC1 changes (~0.1 ppbv) have been de- 
tected at illuminated latitudes. There, SPE-induced active 
chlorine is present in the form of enhanced HOC1 and CIO. 

WACCM4 generally reproduces the SPE-induced variabil- 
ity in the chorine species. Some differences in the geographic 
distribution of these changes arise from the underestima- 
tion of latitudinal redistribution above about 2 hPa. CIO de- 
creases are reproduced quite well at high latitudes, whereas 
the CIO enhancements in the terminator region are underes- 
timated. Moreover, contrary to our expectations, WACCM4 
perturbed and base simulations indicate that this CIO in- 
crease is not due to SPEs. HOC1 variations are modelled very 
well with respect to both magnitude and geographic distribu- 
tion. WACCM4 also reproduces quite well the HC1 depletion 
in the mesosphere (above ~1.5hPa) but underestimates or 
does not show the observed decrease between 1 .5 and 4 hPa. 
CIONO 2 changes observed by MIPAS have also been exam- 
ined. An increase of 0.2 ppbv around 1 hPa in the week after 
the SPE occurrence has been shown. WACCM4 simulations 
with and without proton forcing indicate that this increase is 
dominated by background variability, although SPE-induced 
production might contribute by 0.1 ppbv. 

The most striking differences between satellite data and 
modelling results are (i) the underestimation of the HC1 de- 
pletion mainly below 1 .5 hPa and (ii) the lack of enhanced 
CIO in the terminator region in the WCCM4 simulations. 

The difficulties encountered by Winkler et al. (2011) in 
reproducing the HC1 depletion in stratosphere have been 
confirmed. This model-data discrepancy suggests that defi- 
ciencies in chemical schemes involving chlorine species in 
the upper stratosphere could be present. However the possi- 
ble background variation highlighted in Figs. 4 and 5 could 
also contribute to explain the HC1 discrepancies. Despite 
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the lack of ion chemistry, WACCM4 reproduces (underes- 
timates) mesospheric HC1 changes induced by the SPE of 17 
January at latitudes higher (lower) than 70° N. Since lateral 
mixing is probably underestimated by WACCM4, the real- 
istic reproduction of HC1 concentrations at higher latitudes 
may be the result of compensating errors i.e., too little deple- 
tion, caused by the lack of ion chemistry, plus too little mix- 
ing. Therefore we cannot exclude the importance of the ion 
chemistry in the mesosphere. Addition investigations, taking 
into account the energetic electrons ionization, are necessary. 
Mesospheric HC1 changes induced by the SPE of 20 January 
appear to be slightly overestimated at higher latitudes. Since 
a prior work (Verronen et ah, 2011a) highlighted the impor- 
tance of the mesospheric horizontal mixing on 21 January, 
this is a further element suggesting that lateral mixing could 
be underestimated in WACCM4 during this period. 

The chlorine activation under the SPEs of January 2005 
has been shown. In addition, an HOC1 enhancement accom- 
panied by a similar decrease in HC1 was observed by MLS 
also during the SPEs of December 2006 (i.e., under similar 
seasonal conditions), further corroborating the hypothesis of 
chlorine activation under SPEs. Hints for a deactivation later 
after the SPEs of January 2005 have also been highlighted; 
however, WACCM4 did not fully reproduce this, so no clear 
conclusions can be drawn. More accurate simulations of HO x 
under SPE conditions (Funke et ah, 2011; Jackman et ah, 
2011) could clarify these model-data differences in CIO in 
the the terminator. 

In general, WACCM4 simulations reproduce better the im- 
pact of SPEs on chlorine species for the January 2005 SPE 
than for the 2003 Halloween storm (Funke et ah, 2011). 
The major reason for this difference is that the vortex was 
well established, coincident with the terminator, in January 
2005, while the early winter vortex in October 2003 was par- 
tially illuminated. However, also during January 2005, the 
WACCM4 simulations did not fully capture all the measured 
changes and future work will be necessary to understand and 
remove remaining minor model-data inconsistencies. 
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